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Abstract 

We present the complete analytical result for the two-loop logarithmically enhanced 
contributions to the high energy asymptotic behavior of the vector form factor and the 
four-fermion cross section in a spontaneously broken SU{2) gauge model. On the 
basis of this result we derive the dominant two-loop electroweak corrections to the 
neutral current four-fermion processes at high energies. Previously neglected effects 
of the gauge boson mass difference are included through the next-to-next-to-leading 
logarithmic approximation. 
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1 Introduction 



Recently a new wave of interest in the Sudakov asymptotic regime [1,2] has risen in connec- 
tion with higher-order corrections to electroweak processes at high energies [3-21]. Exper- 
imental and theoretical studies of electroweak interactions have traditionally explored the 
range from very low energies, e.g. through parity violation in atoms, up to energies compara- 
ble to the masses of the W- and Z-bosons at LEP or the Tevatron. The advent of multi-TeV 
colliders like the LHC during the present decade or a future linear electron-positron collider 
will give access to a completely new energy domain. Once the characteristic energies ^/i are 
far larger than the masses of the W- and Z-bosons, My/^zi exclusive reactions like electron- 
positron (or quark-antiquark) annihilation into a pair of fermions or gauge bosons will receive 
virtual corrections enhanced by powers of the large electroweak logarithm ln(s/M^^). 

In Refs. [9,10] we have extended the leading logarithmic (LL) analysis of [8]. The next- 
to-leading logarithmic (NLL) and next-to-next-to-leading logarithmic (NNLL) corrections to 
the high energy asymptotic behavior of the neutral current four-fermion processes have been 
resummed to all orders using the evolution equation approach. Only the light quark case 
was considered and the mass difference between the neutral and charged gauge bosons was 
neglected. On the basis of this result the logarithmically enhanced part of the phenomeno- 
logically important two-loop corrections to the total cross section and to various asymmetries 
was obtained including the ln"(s/M^^) terms with = 2, 3, 4. The results up to NLL 
have been confirmed by the explicit one-loop [5, 12, 13] and two-loop [14-17, 20] calcula- 
tions. The subleading logarithms in the TeV region are comparable to the leading terms due 
to their large numerical coefficients. Thus, the calculation of the remaining two-loop lin- 
ear logarithms is necessary to control the convergence of the logarithmic expansion. These 
corrections represent the next-to-next-to-next-to-leading logarithmic (N^LL) contribution. 
They are of special interest both from the phenomenological and conceptual point of view 
because, in contrast to the higher powers of the logarithm, they are sensitive to the details 
of the gauge boson mass generation. The first results beyond the NNLL approximation have 
been obtained in [16,21,29]. 

In this paper we extend our previous analysis and complete the calculation of the two- 
loop logarithmic corrections to the neutral current four-fermion processes. We incorporate 
previously neglected effects of the gauge boson mass difference. The two-loop logarithmic 
terms are derived within the expansion by regions approach [22-25] by inspecting the struc- 
ture of singularities of the contributions of different regions. The calculation is significantly 
simplified by taking the exponentiation of the logarithmic corrections in the Sudakov limit 
into account. This property naturally appears and can be fully elaborated within the evo- 
lution equation approach [26-28]. To identify the pure QED infrared logarithms which are 
compensated by soft real photon radiation we combine the hard evolution equation which 
governs the dependence of the amplitudes on s with the infrared evolution equation [8] which 
describes the dependence of the amplitude on an infrared regulator. The main results of the 
present paper have been reported earlier in a short letter [29]. 

In Section [2] we present the complete result for the two-loop logarithmic corrections to 
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the form factor which describes fermion scattering in an external Abehan field. Two models 
are considered: a non-Abelian SU{2) gauge theory with the Higgs mechanism of gauge boson 
mass generation and an Abelian U{1) x U{1) theory with two gauge bosons of essentially 
different masses [21]. In Section [3] we generalize the result to the four-fermion process. 
Finally, in Section HI we apply this result to electroweak processes. A brief summary and 
conclusions are given in Section [51 

2 Abelian form factor in the Sudakov limit 

The vector form factor J-' determines the fermion scattering amplitude in an external Abelian 
field. It plays a special role since it is the simplest quantity which includes the complete 
information about the universal coUinear logarithms. This information is directly applicable 
to a process with an arbitrary number of fermions. The form factor is a function of the 
Euclidean momentum transfer = — (pi —P2Y where pi^2 is the incoming/outgoing fermion 
momentum. In the next two sections we consider two characteristic examples: (i) the SU{2) 
gauge model with gauge bosons of a nonzero mass M which emulates the massive gauge 
boson sector of the standard model and (ii) the U{1) x f/(l) gauge model with two gauge 
bosons of essentially different masses which emulates the effect of the Z — 7 mass gap. We 
focus on the asymptotic behavior of the form factor in the Sudakov limit M/Q ^ 1 with 
on-shell massless fermions, Pi = P2 = ^■ 

In the Sudakov limit the coefficients of the perturbative series in the coupling constant 
can be expanded in M^/Q^. To compute the leading term of the series in M^/Q^ we use 
the expansion by regions approach [22-25]. It is based on separating the contributions of 
dynamical modes or regions characteristic for different asymptotic regimes and consists of 
the following steps: 

(i) consider various regions of a loop momentum k and expand, in every region, the inte- 
grand in a Taylor series with respect to the parameters considered small in this region; 

(ii) integrate the expanded integrand over the whole integration domain of the loop mo- 
menta; 

(iii) put to zero any scaleless integral. 

In step (ii) dimensional regularization with d = 4 — 2e space-time dimensions is used to 
handle the divergences. In the Sudakov limit under consideration the following regions are 
relevant [30-32]: 



hard (h) 

1- collinear (Ic) 

2- collinear (2c) 

soft (s) 



k^ r^Q, k+r^ M^/Q , kr^ M, 
kr^ M , 



(1) 
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where k± — ko ± ks, k — (/ci,/c2) and we choose pi,2 = (Q/2, 0, 0, ^(5/2) so that 2pip2 = 
— —s. By k Q, etc. wc mean that every component of k is of order Q. The expansion 
procedure can be facihtated by a technique based on MeUin-Barnes representation which 
can be used not only for evaluating Feynman integrals but also to pick up contributions of 
regions (see Sect. 4.8 of [33] for a general discussion and [34] where the MB representation 
was used for this purpose.) 



2.1 SU{2) model with massive gauge boson 

Let us apply the method to compute the corrections to the form factor in the SU{2) model. 
In one loop the expansion by regions leads to the following decomposition 



^(i)=^W+^(i) + ^(i), 



(2) 



where the subscript c denotes the contribution of both coUinear regions, for a perturbative 
function f{a) we define 

/w=i:(s)"/<"'. (3) 

and the form factor in the Born approximation is normalized to 1. The hard region contri- 
bution, which we will later need, reads 



77(1) _ 



Q 



2e 



(4) 



where Cp = {N^ - l)/(2iV) for a SU{N) gauge group, C(3) = 1.202057 ... is the value of the 
Riemann's zeta-function and all the power-suppressed terms are neglected. For convenience 
we do not include the standard factor (47re~'''^'^(/i^))'' per loop, where 7e — 0.577216. . . is 
Euler's constant. The contributions of all the regions [10] add up to the well known finite 
result 
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where £ = In (Q^/M^). A similar decomposition can be performed in two loops 
The hard-hard part reads 



•J I 



(2) 



hh 
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108 36^ 



2 28 
3^3 + 9^ + 



353 



TT 



/455 7r^\ 1 
^ V108 ^ 36 J e 



Tpris + 



-^+7r^-12C(3) 



TpUf 

+ 0(e') , 



17 

6^3 + 18^2 



Cf 



(5) 



(6) 



(7) 



for a defined in the MS scheme. Here Ca = N,Tf = 1/2, Pq = llCA/'i — 4:TFnf /S — TpUs/S 
is the one-loop beta-function, and Uf (n^) is the number of Dirac fermions (scalars) in 
the fundamental representation. Eq. (I7l) coincides with the known massless two-loop QCD 
result [35, 36] except for the scalar loop contribution which is new. With a renormalized at 
the scale M in the one- loop result, the total two-loop contribution takes the form 





(2) 



Ca + 



(2) 



9^ 9 



Tprif + 



25 
18, 



The coefficients of the second and higher powers of the logarithm are insensitive to the 
infrared structure of the model as explained below. In contrast, the coefficient of the linear 
logarithm does in general depend on the whole mass spectrum of the model. The purely 
Abelian term reads [21] 

-2^ + 2ACi3))cl. (9) 



a(2) 



Massless Dirac fermions give [16] 



a(2) 



-—CfTfUj . 



(10) 



The non-Abelian contribution depends on the details of the gauge boson mass generation. 
For the spontaneously broken gauge group SU{2) with a single Higgs boson of mass Mh = M 
in the fundamental representation explicit calculation leads to 

749 43 o ..... 15 ^ 13 



+ — TT^ - 44C 3 + — VSvr + —V3Ch - 

16 24 ^ ^ 4 ^2 V3 



Here CI2 



1.014942. . . is the value of the Clausen function. For comparison, in the 



(hypothetical) case of a light Higgs boson Mh <^ M this contribution becomes 



A 



(2) 
NA 



A?) 



Zl! + ^.2 _ 44^(3) ^ ^ ^^ci, 

16 48 ^ 8 4 H3 



33 



21 



(12) 



In the electroweak theory inspired model with the SU{2)l gauge group, six left-handed 
fermion doublets, and Mh = M, the result for the two-loop logarithmic corrections reads 



^(2) = 1^4 _ 19^3 

32 48 



+^V3a, 

2 V3 



/463 7 

Vis" 
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61 



15 



, £2 ^ ( 29 - -vr^ - -C(3) + -VStt 



(13) 



The asymptotic dependence of the form factor on Q in the Sudakov limit is governed by the 
hard evolution equation [26-28] 



d 

d\nQ- 



Af2 X 



7(a(x)) + C(a(Q2))+e(a(M2)) 



(14) 
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Its solution is 



^ = Fo(a(M^))exp / — 



X Aj.1 

_7(a(x')) + C(«(x))+e(a(M2)) 



(15) 



By calculating the functions entering the evolution equation order by order in a one gets 
the logarithmic approximations for the form factor. For example, the LL approximation in- 
cludes all the terms of the form and is determined by the one-loop value of 7(a); the 
NLL approximation includes all the terms of the form a"/^^"""^ with m = 0, 1 and requires 
the one-loop values of 7(a), Ci.^^) ^iid ^(a) as well as the one-loop running of a in 7(a); 
and so on. The functions entering the evolution equation can in principle be determined by 
comparing Eq. fll5l) expanded in the coupling constant to the fixed order result for the form 
factor. Within the expansion by regions approach the logarithmic contributions show up as 
singularities of the different regions. One can identify the regions relevant for determining a 
given parameter of the evolution equation and compute them separately up to the required 
accuracy which facilitates the analysis by far. For example, the anomalous dimensions 7(a) 
and C(f^) cire known to be mass- independent and determined by the singularities of the con- 
tribution with all the loop momenta being hard [26-28] . The dimensionally regularized hard 
contribution exponentiates as well [28, 37] with the functions 7(a) and C{oi) parameteriz- 
ing the double and single pole contribution to the exponent. The one- and two-loop hard 
contribution can be written as 



-77(1) 



-r-(2) 
hh 




From Eqs. (|2 H [IID we find 



-2Cf 

Cf 
Cf 



134 



9 

/2545 



2 2\ ^ 40^ 16^ 
-TT I Ca + -^TfUj + —TfUs 

11 . 



I \ TT" 

V 54 9 
'^ + y)TFn.+ g-27r^ + 24C(3))C^ 



(16) 
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At the same time the functions ^(a) and -Fo(«) fix the initial conditions for the evolution 
equation aX Q = M and do depend on the infrared structure of the model. To determine 
the function ^ («) one has to know the singularities of the coUinear region contribution while 
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Fo{a) requires the complete information on the contributions of all the regions. The total 
one- and two-loop form factor can be expressed through the parameters of the evolution 
equation as follows 



-PoC^'^ + Fo^S«)/:^ + (C^^) + + 4'\C^'^ + e^)) c + o{c') . (18) 



With the known values of 7^^'^) and C''^'^'' it is straightforward to obtain the result for the 
remaining functions 

i^'^ = 0, 

c''' = d^Uef + ef + (19) 

(2) 

where the Abehan contribution vanishes =0, the massless Dirac fermions give 



ir={^ + 7y]CFT,n,, (20) 



112 4 

h -> 

27 9 

and for the spontaneously broken SU (2) model with Mh = M we get 

& + e = - 5C(3) + ^73. + ^v^Cl, (0 . (21) 

Note that the functions 7(a) and ^(a) are protected against the Abelian multiloop corrections 
by the properties of the light-cone Wilson loop [27,28,38]. 

The analysis of the evolution equation gives a lot of insight into the structure of the 
logarithmic corrections. For example, Eq. (fTSjl tells us that, up to the NNLL approximation, 
the information on the infrared structure of the model enters through the one-loop coefficients 
^(1) and F^^^ which are insensitive to the details of the mass generation. Thus one can 
compute the coefficient of the two-loop quadratic logarithm with the gauge boson mass 
introduced by hand [10]. 

Note that to complete the N^LL approximation one needs the three-loop value of 7(0;). 
For Us = this has been recently obtained in the context of QCD splitting functions [39] 



7^'^ = 



245 268 2 44 22 A , /836 80 , 112^..A ^ ^ 

(22) 
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+ - 32C(3)) Cj^T^Uf + ^{TFUf)' 
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2.2 U{1) X U{1) model with mass gap 

Let us now discuss the second example, a U{1) x U{1) model with A, a' and M, a for 
masses and coupling constants, respectively. We consider the limit A <^ M and make use of 
the infrared evolution equation which governs the dependence of the form factor JF(A, M, Q) 
on A [8]. The virtual corrections become divergent in the limit A — 0. According to the 
Kinoshita-Lee-Nauenberg theorem [40,41], these divergences are cancelled against the ones 
of the corrections due to the emission of real light gauge bosons of vanishing energy and/or 
collinear to one of the on-shell fermion lines. The singular behavior of the form factor must 
be the same in the full Ua'{l) x ?7a(l) theory and the effective ^/^/(l) model with only the 
light gauge boson. For A -C M -C Q the solution of the infrared evolution equation is given 
by the Abelian part of the exponent (ITS!) with M, a replaced by A, a'. Thus the form factor 
can be written in a factorized form 

^(A, M, Q) = F{M, Q)J'a'i\, Q) + 0{\/M) , (23) 

where JF„/(A,(5) stands for the form factor and F{M,Q) depends both on a and a', 

and incorporates all the logarithms of the form In (Q^/M^). It can be obtained directly by 
calculating the ratio 

"^(A,M,Q)' 



F{M,Q) 



(24) 



Since the function F{M, Q) does not depend on the infrared regularization, the ratio in 
Eq. can be evaluated with A = using dimensional regularization for the infrared 
divergences. The resulting two-parameter perturbative expansion is 



n,m V / 

where 

^(0,0) _ pinfi) ^ ^{0,m) _ jrim) ^ ^36) 

and the two- loop interference term reads [21] 

= (3 - 47r2 + 48C(3)) £ + C(/:°) . (27) 

In the equal mass case, A = M, we have an additional reparameterization symmetry, and 
the form factor is determined by Eq. (fT5l) with the effective coupling a = a' + a so that 
J-'{M, M, Q) = JF^(M, Q). We can now write down the matching relation 

J^(M, M, Q) = C(M, Q)F{M, Q)J'a'{M, Q) , (28) 

where the matching coefficient C{M,Q) represents the effect of the power-suppressed terms 
neglected in Eq. fl23|) . By combining the exphcit results for J^a'{M,Q) and F{M,Q) the 
matching coefficient has found to be C{M,Q) = 1 + 0{a'aC^) [21]. In two-loops it does 
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not contain logarithmic terms, and up to the N^LL accuracy, the product F{M, Q)J^a'{^, Q) 
continuously approaches JF(M, M, Q) as A goes to M. Therefore, to get all the logarithms 
of the heavy gauge boson mass in two-loop approximation for the theory with mass gap, it 
is sufficient to divide the form factor J-'a{M,Q) of the symmetric phase by the form factor 
J-'a'i^, Q) of the effective Ua'{l) theory taken at the symmetric point A = M. Thus we have 
reduced the calculation in the theory with mass gap to the one in the symmetric theory with 
a single mass parameter. The logarithmic terms in the expansion of F{M, Q) exponentiate 
by construction and one can describe the exponent with a set of functions 7(0;, a'), ({a, a'), 
^{a,a') and Fo(a,a') in analogy with Eq. ( |T5l) . The matching procedure can naturally be 
formulated in terms of these functions. For the mass-independent functions we have the all 
order relation 

7(«',") = l{a)-l{a'), 

({a', a) = C(«)-C(«')- (29) 
In two loops we obtain by explicit calculation 

eV,«) = e(«) = o, (30) 

which holds in higher orders for the Abelian model due to the nonrenormalization properties 
discussed in the previous section. Thus, the only nontrivial two-loop matching is for the 
coefficient Fq^'^^ due to the non- logarithmic contribution to C{M,Q) which is beyond the 
accuracy of our analysis. 

Note that the absence of the two-loop linear-logarithmic term in C{M,Q) is an excep- 
tional feature of the Abelian corrections. The general analysis of the evolution equation [10] 
shows that the terms neglected in Eq. (1231) contribute starting from the N^LL approxima- 
tion. Indeed, the solution of the hard evolution equation for F[X, M, Q) which determines 
its dependence on Q is of the form f|T5l) with the infrared sensitive quantities Fq and ^ being 
functions of the ratio A/M. A nontrivial dependence on the mass ratio in general emerges 
first through the two-loop coefficient ^^"^^ due to the interference diagrams with both mas- 
sive and massless gauge bosons. The matching is necessary to take care of the difference 
A/M=i 7^ ^^^^|a/m=o- Thus, for a non- Abelian theory with the mass gap of the standard 
model type, i.e. with interaction between the heavy and light gauge bosons, the matching 
becomes nontrivial already in N^LL approximation. 

3 Four-fermion amplitude 

We consider the four-fermion scattering at fixed angles in the limit where all the kinematical 
invariants are of the same order and are far larger than the gauge boson mass, |s| ~ |t| ~ 
\u\ ^ M^. The analysis of the four-fermion amplitude is complicated by the additional 
kinematical variable and the presence of different isospin and Lorentz structures. We adopt 
the following notation 
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•^LR = i^2Ll^^i^lLi^iRl^^i^3R, (31) 

etc. Here denotes the SU{N) "isospin" generator, pi the momentum of the zth fermion and 
Pi, ps are incoming, and p2, Pa outgoing momenta respectively. Hence t = {pi —p^)'^ = — sx_ 
and u = {pi + PsY = —sx+ where x± = (1 ± cos6')/2 and 9 is the angle between and 
P4. The complete basis consists of four independent chiral amplitudes, each of them of two 
possible isospin structures. For the moment we consider a parity conserving theory, hence 
only two chiral amplitudes are not degenerate. The Born amplitude is given by 

Ab='-^A\ (32) 

s 

The collinear divergences in the hard part of the virtual corrections and the corresponding 
coUinear logarithms are known to factorize. They are responsible, in particular, for the 
double logarithmic contribution and depend only on the properties of the external on-shell 
particles but not on the specific process [26-28,42-45]. This fact is especially clear if a 
physical (Coulomb or axial) gauge is used for the calculation. In this gauge the collinear 
divergences are present only in the self energy insertions to the external particles [28, 43, 
45]. Thus, for each fermion- ant if ermion pair of the four-fermion amplitude the collinear 
logarithms are the same as for the form factor J-' discussed in the previous section. Let us 
denote by A the amplitude with the collinear logarithms factored out. For convenience we 
separate from A all the corrections entering Eq. f[T5]) so that 

^=!^jr2^. (33) 

s 

The resulting amplitude A contains the logarithms of the soft nature corresponding to the 
soft divergences of the hard region contributior0, and the renormalization group logarithms. 
It can be represented as a vector in the color/chiral basis and satisfies the following evolution 
equation [45-47]: 

g^2-^ = xi<ym)A, (34) 

where x(q^) is the matrix of the soft anomalous dimensions. Note that we do not include 
in Eq. fl34l) the pure renormalization group logarithms which can be absorbed by fixing the 
renormalization scale of g in the Born amplitude fl32|) to be Q. The solution of Eq. is 
given by the path-ordered exponent 



A = Pexp 



dx 

—Xia{x)) 

M2 X 



A(a(M2)) ^ (35) 



^ Note that within the expansion by regions the soft x colhnear double logarithmic divergence of the hard 
region contribution is canceled against the ultraviolet divergence of the collinear region i.e. there is no 
one-to-one correspondence between the soft logarithms and the ultraviolet divergences of the soft region 
contributions. 
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where Ao{a) determines the initial conditions for the evolution equation at Q = M. The 
matrix of the soft anomalous dimensions is determined by the coefficients of the single pole 
of the hard region contribution to the exponent ( !35l) . 

In one loop the elements of the matrix x(q^) do not depend on chirality and read [9] 



Xxd 



4 



-2Ca (In (x 

CpTp 
N 



zvr) + 4 Cf 



Tf 
N 



In 



In 



Xdx 

Xdd 



41n^ 



(36) 



In terms of the functions introduced above the one-loop correction reads 



(1) 



where 



^0 



^0 LL-^LL 



A, 



(1)A 



LR-^LR 



(37) 



(3^ 



For the present two-loop analysis of the annihilation cross section only the real part of the 
coefficients Aq'^ is needed (see e.g. [10]), 



Re 
Re 
Re 
Re 
where 



A 



A 



(1)A 
LL 

(1)A 
LR 



85 



— ) /(x+, x.) + Ca[ — + t( ) - —TpUf - -Tpn 
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20, 



Ca 



^0 LL 



CpTp 



A 



LR 



N 

CpTp 



N 



N 

/(x+, x_) , 



^--^\f{x.,X^)+CA{- 



85 
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7r 



20^ 8^ 



f{x+,x. 



2 , X- 
— In a;_ H 



In X. 



X, 



(39) 



(40) 



The two-loop matrix x*^^^ can be extracted from the result for the single pole part of the 
hard contribution to the four-quark amplitude [48-50]. For vanishing beta- function it reads 



(41) 
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After averaging over the isospin Eq. (HTj) agrees with the result of Ref. [51]. An alternative 
approach to compute the soft anomalous dimension matrix is based on studying the renor- 
malization properties of the product of light-like Wilson lines (see, e.g. [52]). This approach 
has been used in Ref. [53] for the calculation of the soft anomalous dimension matrix in two 
loops. The result of Ref. [53] is in agreement with Eq. (14T1) . 

The part of the matrix x^"^^ proportional to the beta-function originates from the running 
of the coupling constant in the one-loop nonlogarithmic term. We find that Eq. fl39|) gets 
contributions only from the hard region so that the effective renormalization scale of a there 
is Q. Hence 



R /'4(1)A \\ I 4(1)A ylA , ~AiX)d Ad , A{l)d A, 

'PO LL-^LL "T ^0 LR-^LR + ^0 LL-^LL + ^0 LR-^LR 



(42) 



i.e. in this order the soft anomalous dimension matrix does depend on the chirality. The 
matrix x*'^'^'* as well as the coefficients ^o^^/j are determined by the hard region contribution 
and are insensitive to the details of the gauge boson mass generation. 

Note that all the nontrivial soft region contributions in the expansion of the multiscale 
two- loop integrals in the Sudakov limit are power-suppressed [23-25]. The unsuppressed 
soft region contributions originate from one-scale self-energy insertions into an external on- 
shell fermion or into a collinear gluon propagator. On the other hand one power of the 
logarithm in the double logarithmic contribution is of soft nature. Thus the two-loop soft 
logarithms in A can be related to the collinear divergences though the latter factorize and 
do not contribute to the soft anomalous dimension matrix themselves. This, in particular, 
explains the presence of the 7''^'* contribution to x^'^'^ ■ 

The two-loop correction to the four-fermion amplitude is obtained by the direct general- 
ization of the form factor analysis with the only complication related to the matrix structure 
of Eq. ([35D 



(1)£3 



A' + (2Fo«x« + xS) A' + (2C« + X«) C + 0{C 



A' 



(43) 



X,dd 



and the matrix structure of the product x'^^^'Aq 



(1) 



where we used the fact that ^^^^ 
is implied. 

Let us again discuss the standard model inspired example considered in the previous 
section. With the result for the amplitudes it is straightforward to compute the one- and 
two-loop corrections to the total cross section of the four-fermion annihilation process. For 
the annihilation process /'/' ff one has to make the analytical continuation of the above 
result to the Minkowskian region of negative = —s according to the s + iO prescription. 
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The above approximation is formally not valid for the small angle region 9 < M/ ^/s, which, 
however, gives only a power-suppressed contribution to the total cross section. For the 
SU{2)l model we obtain 



(2) 



and 



(2) 



lc\s) - ^^C\s) + 



4855 37 



+ 157371 + 2673012 i-]]C 



+ ^n' C'is) + 



48049 1679 



216 



1^ 



vr^ - 122C(3) 



(4.50 C\s) - 74.83 C^{s) + 391.45 C'^{s) - 717.49 C{s)) as 



(44) 



-cHs) - —cH 



125 



799 37 
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-15V37r + 26^3012 f|] ]C{s) 



38005 383 



216 



1^ 



-TT 



122C(3) 



(4.50 C^{s) - 20.83 C^{s) + 32.95 C'^{s) - 53.38 C{s)) ae 



(45) 



for the initial and final state fermions of the same or opposite eigenvalue of the isospin 
operator t^, respectively. Here 0"^ is the Born cross section with the MS coupling constants 
renormalized at the scale ^/s and £(s) = ln(s/M^). 



4 Two-loop electroweak logarithms 

In Born approximation, the amplitude for the neutral current process /'/' ff reads 

-4b = ^ E (r;.r; + 4^)^f:/. (46) 

where 

M'J = I'nMji.fj , (47) 

tw = tanOw with 9\y being the Weinberg angle and Tf (Yj) is the isospin (hypercharge) of 
the fermion which depends on the fermion chirality. The one-loop electroweak correction to 
the amplitude is well known (see [54] for the most general result). The calculation of the 
two-loop electroweak corrections even in the high energy limit is a challenging theoretical 
problem at the limit of available computational techniques. It is complicated in particular 
by the presence of the mass gap and mixing in the gauge sector. In Sect. 14.11 we develop 
the approach of Ref. [21] and reduce the analysis of the dominant two-loop logarithmic 
electroweak corrections to a problem with a single mass parameter which has been solved in 
the previous section. In Sect. 14.21 we present the explicit expression for the two-loop virtual 
linear-logarithmic electroweak correction to the four-fermion amplitudes. In Sect. 14.31 we 
complete the analysis of Refs. [9, 10] by computing the leading effects of the W — Z mass 
splitting through NNLL approximation. In Sect. l4.4l we introduce infrared safe semi-inclusive 
cross sections and in Sect. 14.51 we give numerical results for the two-loop corrections to the 
cross sections and various asymmetries. 
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4.1 Separating QED infrared logarithms 



The main difference between the analysis of the electroweak standard model with the spon- 
taneously broken SUl{'2) x f/(l) gauge group and the treatment of the pure SUl{2) case 
considered above is the presence of the massless photon which results in infrared divergences 
of fully exclusive cross sections. We regularize these divergences by giving the photon a 
small mass A. The dependence of the virtual corrections on A in the limit A^ <^ <^ 
is governed by the QED infrared evolution equation and is given by the factor 



U = f/o(ae(Q')) 



X exp 



47r 



Q} + Qf-^C^N,{Q} + Q},) + fN,ln 





Qf'Qf 



+ 4 



ae 160 



Qf'Qf 



In 



'Q^ 
A2 



vr 27 



{q} + Q},) In^ 




(48) 



where ae is the MS QED coupling constant and Qf is the electric charge of the fermion. The 
NNLL approximation for U can be found in [10] and to derive the two-loop linear logarithm 
in the exponent (HHj) we use the following expressions for the QED parameters: 



iV,+ Q-27r2 + 24C(3))Q 



272 

640 , fx+\ ^ ^ 



(49) 



The expressions for ^^^^ and x^e ^ can be obtained by substituting rij — > 8A^g/3 into the 
general formulae. The coefficient Uo{ae{Q'^)) in Eq. fHSl) is a two-component vector in the 
chiral basis. We have a freedom in definition of this quantity because it does not depend 
on A and can be absorbed into the electroweak part of the corrections. It is convenient 
to normalize the QED factor to U\q=x = 1 so that f/g^^ = 0. Note that the structure of 



contribution to the exponent fHHj) and the NgQ'j part of the linear 



Eq. fHS]) differs from the one of the solution of the hard evolution equations (IT^ [25]) because 
cte in the preexponential factor Uo{ae) is renormalized at Q so that all the dependence on 
A is contained in the exponent. By contrast, in Fo{a) and Ao{a) the coupling constant is 
renormalized at M so the that all the dependence on Q is contained in the exponent. This 
is why there is no X^e^lpo 

logarithm coefficient there differs from the one of (^^^ + .^p-* . 

Our goal now is to separate the above infrared divergent QED contribution from the 
total two-loop corrections to get the pure electroweak logarithms ln(Q^/M^^). Within the 
evolution equation approach [10] it has been found that the electroweak and QED logarithms 
up to the NNLL approximation can be disentangled by means of the following two-step 
procedure: 
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(i) the corrections are evaluated using the fields of unbroken symmetry phase with all the 
gauge bosons of the same mass M ^ Mz^w, i-^- without mass gap; 

(ii) the QED contribution (l48l) with A = M is factorized leaving the pure electroweak 



This reduces the calculation of the two-loop electroweak logarithms up to the quadratic 
term to a problem with a single mass parameter. Then the effect of the Z — W boson 
mass splitting can systematically be taken into account within an expansion around the 
equal mass approximation [21]. In general the above two-step procedure is not vahd in the 
N^LL approximation which is sensitive to fine details of the gauge boson mass generation. 
For the exact calculation of the coefficient of the two-loop linear-logarithmic term one has 
to use the true mass eigenstates of the standard model. The evaluation of the corrections 
in this case becomes a very complicated multiscale problem. The analysis, however, is 
drastically simplified in a model with a Higgs boson of zero hypercharge. In this model the 
mixing is absent and the above two-step procedure can be applied to disentangle all the two- 
loop logarithms of the S't/i(2) gauge boson mass from the infrared logarithms associated 
with the massless hypercharge gauge boson (see the discussion below). In the standard 
model the mixing of the gauge bosons results in a linear-logarithmic contribution which 
is not accounted for in this approximation. It is, however, suppressed by the small factor 
sin^ 6w = s^r ~ 0.2. Therefore, the approximation gives an estimate of the coefficient in 
front of the linear electroweak logarithm with 20% accuracy. As we will see, a 20% error in 
the coefficient in front of the two-loop linear electroweak logarithm leads to an uncertainty 
comparable to the nonlogarithmic contribution and is practically negligible. If we also neglect 
the difference between and Mz^w^ the calculation involves a single mass parameter at 
every step and the results of the previous sections can directly be applied to the isospin 
SU{2)l gauge group with the coupling g and the hypercharge U{1) gauge group with the 
coupling twQ- 

The procedure of separating the QED contribution can naturally be formulated in terms 
of the functions parameterizing the solution of the hard evolution equation for the four- 
fermion amplitude. The functions 7(a), C(q^)) and xi'^) are mass- independent. Therefore, 
these functions parameterizing the electroweak logarithms can be obtained by subtracting 
the QED quantities from the result for the unbroken symmetry phase to all orders in the 
coupling constants (c/. Eq. fl29l) ) without referring to the simplified model. To the order of 
interest they can be found in [10] or easily derived from the results of the previous sections. 
For example, after the subtraction, the /^o = part of the matrix x*'^^ for I, J = L takes the 



logarithms. 



form 




/3o=0 



( 




130 27r2\ 



9 3 J 



4 (In (x+) +i7!-) + 2 In 
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Xxd 



(2) 



(2) 
/(.dd 



/3o=0 



/3o=0 



/3o=0 



80,, 520 87r2\, , „^ 
-—Na + - ^ 1 In I 

X- 




NgSwQf'Qf 



In 



For / or/and J = i? it is reduced to 



/3o=0 




(50) 



(51) 



Note that no QED subtraction is necessary for x \po because of the specific normahzation 
of the QED factor U. Thus ^^"^^ is the only two- loop coefficient at the order of interest which 
is sensitive to fine details of the gauge boson mass generation. We evaluate it approximately 
by using the above simplified model with the Higgs boson of zero hypercharge and of the 
mass Mh = Mw,z- In this model the interference diagrams including the heavy SUii^) and 
the light hypercharge U{1) gauge bosons are identical to the ones of the U{1) x f/(l) model 
discussed in Sect. 12.21 and the corresponding contribution to vanishes. At the same time 
the pure non-Abelian contribution to is given by Eq. f|T9l) . Collecting all the pieces for 
the two-loop contribution to the scheme-independent sum ({a)+C,{<y) in the simplified model 
we obtain 



^(2) + ^(2) 



34,, 
3 ^ 



749 43 



176 



12 



-TC 



C(3) 



5v^vr + |y3Cl2(|) 



+ (3-47r2 + 48C(3)) 4^ 



(52) 



where Tj = and 1/2 for the right- and left-handed fermions, respectively. 



4.2 Two-loop virtual corrections 

The expressions for the two-loop corrections to the electroweak amplitudes is obtained by 
projecting the result of Sect. [3] on a relevant initial/final state with the proper assignment 
of isospin/hypercharge. For example, the projection of the basis (|3T1) on the states corre- 
sponding to the neutral current processes reads A^j TfTf,A{'/, Ajj A{'/ . The only 
complication in combinatorics is related to the fact that now we are having different gauge 
groups for the fermions of different chirality. To take into account the light fermion contri- 
bution one has to replace nj ^ 2Ng for the SU{2)l and nf 5Ng/3 for the hypercharge 
f/(l) gauge group, where A^^ = 3 stands for the number of generations. The result for the 
linear-logarithmic corrections is given in the following approximation: 

(i) the fermion coupling to the Higgs boson is neglected, i. e. all the fermions are considered 
to be massless; 
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(ii) the W — Z mass splitting is neglected; 

(iii) the approximation (|52l) is used for the two- loop coefficient C*-^-* + 

The accuracy of the approximation is discussed in Sect. (14. 5p . The result for the n-loop 
correction to the amphtude ( H6l) can be decomposed as 

j^H ^ _|_ _)_ a\^]sill + • • • • (53) 

Explicit expressions for the corrections up to A^^^j^^ can be found in [9,10]. The two-loop 
N^LL term is new. For convenience we split it in parts as follows: 



N^LL 



^I'^N^LL ■ 



(54) 



i=l 



The correction corresponding to the C^^"* ^^id Fg^''^'^^^ terms of Eq. fH5]) is 



a^A.A^XL 



+^V3Ch f ^) + (3 - 47r2 + 48C(3)) 4^ ) x T;(T^ + 1) + (/ /') 



E 
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X 



T^/(T^/ + 1)+ 4^ + (/--/') 



l{Q')aAg 



(55) 



where a = (jf^/lGvr^ and /(Q^) = aln (Q^/M^). The correction corresponding to the x^^"* term 
of Eq. (H3l) for vanishing beta-function reads 



A2A 
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(56) 



while the /3o part of x^"^^ term gives 
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(57) 



where 5ij = 1 for I = J and zero otherwise. Note that our analysis imphes the MS couphng 
constants renormahzed at the scale M in the one-loop result with the exception of the one 
power of the coupling constants originating from the Born amplitude (l46ll renormalized at 
the scale Q. The correction corresponding to the C^^-'Re 
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The correction corresponding to the x^^^^o^^ term of Eq. (H3!) reads 

-4 (In (x+) + + In (^^^ (2 + ^F/^F/) j T;,T/ 

+ In - 4.^^g,,Q,) (t3t3 + | 

/(Q2)a^f7. (60) 
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The above expressions are derived for Euclidean positive . For the annihilation process one 
has to make the analytical continuation to the Minkowskian region of negative = — (s+iO) 



and it is natural to normalize the QED factor to L{\ 



=M2 



1 at the Minkowskian point 



rather than at Q'^ 



M^. If we use this normalization condition then Re 
-7r^((5^ + Qj/) and no QED subtraction is necessary for (the real part of) the amp 



itude at 



4.3 The effect ofW — Z mass splitting 

The NLL and NNLL result for the electroweak corrections has been obtained in Refs. [9, 10] 
neglecting the Z — W boson mass splitting which is suppressed by a small factor 5m = 
1 — My^/Ml 0.2. The effect of the mass splitting in the NLL and NNLL terms can, 
however, be comparable to the N'^LL contribution and should be taken into account at this 
accuracy level. In the NLL and NNLL approximation the mass splitting can easily be taken 
into account through the modification of the one- loop parameters Fq'\ and JiQ^ Let us 
take Mw as the argument of the logarithms. Then the mass splitting corrections originate 
from the Z boson contribution. From the explicit result for the one-loop diagrams in the 
leading order in ancl^ 5m oc we obtain 
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By using Eq. fl6T|) it is straightforward to get the leading effect of the mass splitting in one 
and two loops through the NNLL approximation. The corrections to the amplitudes read 
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We do not express 5m through to emphasize the origin of the corrections. 
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Yf'Yf 
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(62) 



which completes the result of Refs. [9, 10]. Note that the exponentiation of the mass split- 
ting contribution due to ^.^'■^^ has been observed first by the explicit calculation in NLL 
approximation [20]. 



4.4 Semi- inclusive cross sections 

To get the infrared safe result one has to take into account the real photon emission in 
an inclusive way. In practice, the massive gauge bosons are supposed to be detected as 
separate particles. Thus it is of little physical sense to treat the hard photons with energies 
far larger than Mz^w separately because of gauge symmetry restoration. In particular, the 
radiation of the hard real photons is not of the Poisson type because of its non-Abelian 
SU{2)l component. Therefore, we restrict the analysis to semi-inclusive cross sections with 
the real emission only of photons with energies far smaller than Mz,w, which is of pure 
QED nature. To derive the result for such a cross section one has to add to the expressions 
given above the QED corrections due to the real photon emission and the pure QED virtual 
corrections which are determined for mj <^ \ <^ M, where is a light fermion mass, by 
Eqs. fl48p . However, in standard QED applications a nonzero mass mj is kept for an on-shell 
fermion which regularizes the collinear divergences. To derive the QED factor for A far less 
than the fermion mass X <^ rrif <^ M one has to change the kernel of the infrared evolution 
equation and match the new solution to Eq. fHHj) at the point X = rrif. For several light 
flavors of significantly different masses the matching is necessary when A crosses the value 
of each fermion mass. The sum of the virtual and real QED corrections to the cross section 
gives a factor which depends on s, the fermion masses, and the experimental cuts, but not 
on Mz,w- The detailed analysis of the QED corrections goes beyond the scope of the present 
paper. In the single flavor case the second order QED corrections including the soft real 
emission are known beyond the logarithmic approximation (see [55] and references therein). 



4.5 Numerical estimates 

With the expressions for the chiral amplitudes at hand, we can compute the total two-loop 
logarithmic virtual corrections to the basic observables for the neutral current four-fermion 
processes. Let us consider the total cross sections of the quark-antiquark/ /i"*"/!" production 
in e~^e~ annihilation. In one and two loops the logarithmic corrections read 

R{e^e--^QQ) = 1 - 1.66 L(s) + 5.60 /(s) - 8.39 a 

+ 1.93 L2(s) - 11.28 L(s)/(s) + 33.79 /2(s) - 60.87/(s) a, 
R{e^e-^qq) = 1 - 2. 18L(s) + 20.94/(s) - 35.07a 

+ 2.79 L^{s) - 51.98 L{s) /(s) + 321.20 /^(s) - 757.35 l{s) a , 
i?(e+e- ^ /x+/i") = 1 - 1.39 L{s) + 10.35 l{s) - 21.26 a 

+ 1.42 L2(s) - 20.33 L(s)/(s) + 112.57 /^(s) - 312.90 /(s) a , (63) 
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where L(s) = aln^(s/M^), Q = u,c,t, q = d,s,b, R{e^e — > jj^fi ) = a/aB{e^e 
/i"*"/!") and so on. The MS couphngs in the Born cross section are renormahzed at ^/s. 
Numerically, we have L{s) = 0.07 (0.11) and l{s) = 0.014 (0.017) for = 1 TeV and 2 TeV, 
respectively. Here M = Mw has been chosen for the infrared cutoff and a = 2.69 ■ 10~^, 
Syy = 0.231 for the MS couplings renormalized at the gauge boson mass. The complete 
one-loop corrections are known exactly and we have included the dominant one-loop terms 
in Eqs. fl63l) - fl67|) to demonstrate the structure of the expansion rather than for precise 
numerical estimates though the above expressions approximate the exact one-loop result 
with 1% accuracy in the TeV region. In Eq. (1551) we included the leading correction in the 
W — Z mass difference through NNLL approximation. The coefficient of the linear logarithm 
is computed in the approximation of Sect. (14. ip . In the case of a quark-antiquark final state 
the strong interaction could also produce logarithmically growing terms. For massless quarks 
the complete 0{al) corrections including the bremsstrahlung effects can be found in [35]. 
Note that Eqs. (I63p are symmetric under exchange of the initial and final state fermions 
and, therefore, also applicable to the Drell-Yan processes at hadron colliders. In this case 
the two-loop QCD corrections are given in [36]. For the total cross sections of the four-quark 
electroweak processes we obtain 

R{Q'Q' QQ) = 1 - 2.07L(s) + 19.03 /(s) - 32.63 a 

+ 2.67L2(s) -46.64L(s)/(s) + 278.94 /2(s) - 666.05 /(s) a, 
R{QQ^qq) = 1 -2.56 L(s) + 8.49 /(s) - 11.94 a 

+ 3.53 L^(s) - 20.39 L{s) l{s) + 65.20 l^{s) - 91.92 l{s) a , 
Riq'q'-^qq) = 1 - 2.87 L(s) + 25.63 /(s) - 38.89 a 

+ 4.20 L2(s) -71.87L(s)/(s) + 423.61 /^(s) -919.35 /(s) a. (64) 

Our results can easily be generalized to // f f processes with identical quarks by including 
the t-channel contribution. 

For e~^e~ annihilation we also give a numerical estimate of corrections to the cross section 
asymmetries. In the case of the forward-backward asymmetry A^^ (the difference of the cross 
section integrated over forward and backward hemispheres with respect to the electron beam 
direction divided by the total cross section) we get 

l-0.09L(s)-1.22/(s) + 1.77a 
+ 0.12 L^{s) + 0.59 L{s) l{s) - 1.65 P{s) + 3.39 /(s) a , 
1 - 0.14L(s) + 7.17/(s) - 10.84a 

+ 0.02L2(s) - 1.32L(s)/(s) - 33.07/2(s) + 93.44/(s)a, 
1 - 0.04 L{s) + 5.50 l{s) - 14.43 a 

+ 0.27L2(s) -6.43 L(s)/(s) + 22.91 /2(s) -33.11 /(s) a, (65) 

where R^^ = A^^ /A^^. For the left-right asymmetry A^^ (the difference of the cross 
sections of the production of left- and right-handed particles divided by the total cross 



/?^^(e+e- ^ QQ) = 
R^^{e+e-^qq) = 
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(a) 



(b) 



Figure 1: Separate logarithmic contributions to R{e^e — > qq) in % to the Born approximation: 
(a) the one-loop LL (In^is/M^), long-dashed line), NLL {ln^{s/M'^), dot-dashed line) and N^LL 
(ln°(s/M2), solid line) terms; (b) the two-loop LL {ln^{s/M'^), short-dashed line), NLL {ln^{s/M'^), 
long-dashed line), NNLL (ln^(s/M^), dot-dashed line) and N^LL (ln^(s/M^), solid line) terms. 



section) we obtain in the same notations 



e ' e 



R 



LR, 



e e 



QQ) = 1 -4.48 L(s) + 17.51 /(s) - 13.16 a 

- 1.16 L^{s) + 15.66 L{s) l{s) - 43.50 /^(s) + 44.05 l{s) a , 
^qq) = 1 - 1.12L(s) + 12.05/(s) - 16.44a 

- 0.81 L2(s) + 18.02 L{s) l{s) - 130.74 f{s) + 278.71 l{s) a , 
ix^ix') = 1 - 13.24L(s) + 116.58/(s) - 148.42a 

- 0.79 L2(s) + 23.68 L(s) l{s) - 155.46/^(5) - 116.67/(s) a. 



(66) 

Finally, for the left-right asymmetry A^^ (the difference of the cross sections for the left- 
and right-handed initial state particles divided by the total cross section) which differs from 
j^LR. quark- ant iquark final state we have 



R 



LR, 



e e 



QQ) 



R 



LR, 



e e 



- 1 -2.75L(s) + 10.60/(s) -9.05a 

- 0.91 L2(s) + 11.16 L{s) l{s) - 33.49 l^{s) + 28.28 /(s) a , 
qq) = 1 - 1.07L(s) + 11. 75/(s) - 16.21a 

- 0.77 L^{s) + 17.06 L(s)/(s) - 125.18 /^(s) + 267.60 /(s) a . 



(67) 

The numerical structure of the corrections in the case of e+e^ annihilation is shown in 
Figs. [Ilini In Fig. [1] the values of different logarithmic contributions to R{e^e~ — > qq) are 
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(a) 




(b) 



Figure 2: The total logarithmic corrections to -R(e+e~ QQ) (dashed line), R{e'^e~ qq) (dot- 
dashed line) and R{e^e~ IJ'~^IJ'~) (solid line) in % to the Born approximation: (a) the one-loop 
correction up to N^LL term; (b) the two-loop correction up to N^LL term. 

plotted separately as functions of s. In the 1 — 2 TeV region the two-loop LL, NLL, NNLL 
and N^LL corrections to the cross sections can be as large as 1 — 4%, 5 — 10%, 5 — 10% 
and 2 — 3% respectively. The two-loop logarithmic terms have a sign-alternating structure 
resulting in significant cancellations typical for the Sudakov limit [10,16,21]. Although the 
individual logarithmic contributions can be as large as 10%, their sum does not exceed 1% 
in absolute value at energies below 2 TeV for all the cross sections (see Fig. [2]). In the 
region of a few TeV the corrections do not reach the double-logarithmic asymptotics. The 
quartic, cubic and quadratic logarithms are comparable in magnitude. Then the logarithmic 
expansion starts to converge. Still, the linear-logarithmic contribution must be included to 
reduce the theoretical uncertainty below 1%. The two-loop logarithmic corrections to the 
asymmetries also amount up to 1% in absolute value at energies below 2 TeV with the only 
exception of R^^{e~^e^ l^^l^ ) (see Fig. [3]). 

Let us discuss the accuracy of our result. On the basis of the explicit evaluation of the 
light fermion/scalar [16] and the Abelian contribution [21] we estimate the uncalculated two- 
loop nonlogarithmic term to be at a few permill level. The power-suppressed terms do not 
exceed a permill in magnitude for ^/s > 500 GeV as well [16]. By comparing our numerical 
estimates to the equal mass approximation given in the Appendix we find that the leading 
effect of the W — Z mass splitting results in a variation of at most 5% of the coefficients of the 
two-loop cubic and quadratic logarithms. Thus the expansion in the W — Z mass difference 
converges well for these coefficients and the leading correction term taken into account in our 
evaluation is sufficient for a permill accuracy of the cross sections. The main uncertainty of 
our result for the two-loop single logarithmic contribution is due to our approximation for the 

coefficient. Neglecting the gauge boson mixing effects, which are suppressed by a factor of 
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(a) 



(b) 



Figure 3: (a) The total logarithmic two- loop corrections to the forward-backward asymmetry 

R^^{e+e- QQ) (dashed line), R^^{e+e- qq) (dot-dashed line) and R^^{e+e- 

(solid line) in % to the Born approximation, (b) The same as (a) but for the left-right asymmetry 

sin^ 9\Y, brings an error of 20% to the coefficient of the two- loop single logarithm. Neglecting 
the difference between the Higgs and gauge boson masses does not lead to a numerically 
important error because the scalar boson contribution is relatively small. By comparing the 
results for Mh = Mw,z and Mh -C Mw,z we estimate the corresponding uncertainty in the 
coefficient of the two- loop single logarithm to be about 5%. The same estimate is true for the 
uncertainty due to the top quark mass effect on the ti virtual pair contribution. Hence for the 
production of light fermions our formulae are supposed to approximate the exact coefficients 
of the two-loop linear logarithms with approximately 20% accuracy which results in a few 
permill uncertainty in the cross sections. By adding up the errors from different sources 
in quadrature we find the total uncertainty of the cross section to be in a few permill - 
one percent range depending on the process. Thus we get an accurate estimate of the two- 
loop correction, which is sufficient for practical applications to the future collider physics. 
The only essential deviation of the exact two-loop logarithmic contributions from our result 
Eqs. (!63l - [67l) for the production of the third generation quarks is due to the large top quark 
Yukawa coupling. The corresponding corrections are known to NLL approximation and can 
numerically be as important as the generic non- Yukawa ones [11,13,56,57]. 

5 Summary 

In the present paper we have derived the analytical result for the two-loop logarithmic cor- 
rections to the vector form factor and the four-fermion cross sections in the spontaneously 
broken SUl{2) model by combining the explicit calculation and the evolution equation ap- 
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proach. We have completed the analysis of the dominant logarithmically enhanced two-loop 
electroweak corrections to the basic observables of the neutral current four-fermion processes 
at high energy. The two-loop linear-logarithmic contribution has been obtained with an es- 
timated accuracy of 20%. The W — Z mass sphtting effect neglected in [9,10] has been taken 
into account through the NNLL approximation. Our result for the two-loop logarithmic 
corrections along with the known exact one-loop expressions approximates the cross sections 
with an accuracy of a few permill, which is sufficient for practical applications to the future 
collider physics. The general approach for the calculation of the logarithmically enhanced 
corrections developed in the present paper can also be applied to the analysis of the gauge 
boson production and the supersymmetric extensions of the standard model where only the 
NLL approximation is available so far [18, 19]. 
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Appendix 

The NNLL result for the cross sections and asymmetries in the approximation Mz = My/ 
can be obtained from the expressions of Ref. [lojf] 

1 - 1.66 L(s) + 5.31 /(s) -8.36 a 
+ 1.93 L^{s) - 10.59 L{s)l{s) + 31.40 /^(s) , 
1 - 2.18 L{s) + 20.58 l{s) - 34.02 a 
+ 2.79 L2(s) - 51.04 L{s)l{s) + 309.34 t^^s) , 
1 - 1.39 L{s) + 10.12 l{s) - 20.61 a 
+ 1.42 L2(s) - 19.81 L{s)l{s) + 107.03 /^(s) , (68) 

1 - 0.09 L(s) - 1.23 /(s) + 1.47a 
+ 0.12L2(s) + 0.64L(s)/(s) - 1.40 /^(s) , 
1 - 0.14L(s) + 7.15/(s) - 10.43a 
+ 0.02 L2(s) - 1.31 L{s)l{s) - 33.46 /^(s) , 
1 - 0.04 L{s) + 5.49 l{s) - 14.03 a 
+ 0.27L2(s) -6.32L(s)/(s) + 21.01 /^(s), (69) 



R{e+e-^QQ) = 
R{e^e~ qq) = 



i?^^(e+e-^gg) = 

R^^{e+e- qq) = 
i?^^(e+e- ^ /i+/i-) = 



R^^{e+e- ^ QQ) = 1 - 4.48 L(s) + 16.66 /(s) - 13.28 a 

- 1.16 L2(s) + 15.21 L{s)l{s) - 41.79 f{s) , 
R^^ie^e- ^qq) = 1-1.12 L{s) + 11.86 l{s) - 15.83 a 

- 0.81 L\s) + 17.74 L{s)l{s) - 127.05 P{s) , 
i?^-^(e+e- ^ /i+/i") = 1 - 13.24L(s) + 113.77/(s) - 139.94a 

- 0.79 L'^is) + 23.34 L{s)l{s) - 155.36 /^(s) . (70) 



■^Throughout Sect 4. of Ref. [10] the terms with the factor {aNg + b)t^f, where a and b stand for some 
constants, should be multiphed by an extra t^r, and the terms with the factor Ngs'^^ should be multiplied 
by an extra s^. This results in a small change of the numerical estimates. 
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R^^(e+e- QQ) = 1 - 2.75 L(s) + 10.07 - 9.02 a 

- 0.91 L'^(s) + 10.80 L(s)l(s) - 32.10 f(s) , 
R^^{e+e- ^qq) = 1 - 1.07 L{s) + 11.56 - 15.60 a 

- 0.77 L^(s) + 16.78 L(s)l(s) - 121.56 f(s) . (71) 
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